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Molecular  dynamics  is  used  to  model  the  structure  and  dynamics  of  electric  double  layers 
at  a  charged  metal  surface  and  to  ask  questions  of  spectroscopic  interest.  In  particular: 
What  is  the  molecular  basis  for  changing  the  position  of  the  outer  Helmholtz  plane 
(OIIP)  thereby  tuning  the  electric  Held  across  the  inner  layer?  Another  topic  is  the  mi¬ 
croscopic  basis  of  modulation  spectroscopies  like  SNIP  MRS  Out  Ml)  simulations  of  an 
immersed  electrode  show  features  corresponding  to:  compact  layer,  diffuse  layer,  highly 
oriented  water  layer  next  to  the  metal  when  the  electrode  is  charged  and  ions  are  present, 
penetration  of  nominally  diffuse  layer  species  into  inner  layer,  ion  pair  formation  between 
contact  adsorbed  ion  and  diffuse  layer  ion  when  the  electrode  is  uncharged,  poorly  ori¬ 
ented  surface  water  when  the  electrode  is  uncharged.  All  these  properties  arise  from  the 
model  with  the  restriction  that  charge  on  the  metal  and  aqueous  phase  sums  to  zero,  ic., 
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Introduction  and  Model 

Wc  use  molecular  dynamics  (MD)  to  model  the  structure  and  dy¬ 
namics  of  electric  double  layers  at  a  charged  metal  surface  and  to 
ask  questions  of  spectroscopic  interest.  In  particular:  What  is  the 
molecular  basis  for  changing  the  position  of  the  outer  Helmholtz, 
plane  (OIIP)  thereby  tuning  the  electric  field  *•  2  across  the  inner 
layer?  Another  topic  is  the  microscopic  basis  of  modulation 
spectroscopies  like  SNIFTIRS  3.  Our  MD  simulations  of  an  im¬ 
mersed  electrode  show  features  corresponding  to:  compact  layer, 
diffuse  layer,  highly  oriented  water  layer  next  to  the  metal  when 
the  electrode  is  charged  and  ions  are  present,  penetration  of  nom¬ 
inally  diffuse  layer  species  into  inner  layer,  ion  pair  formation  be¬ 
tween  contact  adsorbed  ion  and  diffuse  layer  ion  when  the 
electrode  is  uncharged,  poorly  oriented  surface  water  when  the 
electrode  is  uncharged.  All  these  properties  arise  from  the  model 
with  the  restriction  that  charge  on  the  metal  and  aqueous  phase 
sums  to  zero,  ic.,  qM  Uqions=0. 

Interaction  with  the  metal  was  represented  by  a  9-3  potcntal  for 
Pauli  repulsion  and  attractive  dispersive  interactions,  and  an  im¬ 
age  potential  for  interaction  with  the  conduction  electrons.  On 
the  side  opposite  the  metal  the  electrolyte  was  constrained  by  the 
9-3  potential  of  a  'dielectric-like’  bounding  surface.  This  boundary 
limited  the  extent  of  the  fluid  phase,  and  kept  the  calculations 
tractable.  System  composition  was  mM  1  +nX M  (N-m-n)IIjO 


where  N  is  the  number  of  water  molecules  in  the  absence  of  ions, 
and  (m,n)  =  (0 ,()),( I,  0),(0,  1),(1,  I ),( 2,  I) .  I  Icinzingcr  parame¬ 

ters  4  for  st2  water  mode!  and  alkali  halides  MX  were  used. 
Electrostatic  fields  were  calculated  exactly  (no  cut  offs)  by  the  fast 
multipolc  method. 

Double  Layer  on  Charged  Electrode 

We  display  a  sample  result  here  (see  Figure  I)  for  a  system  com¬ 
posed  of  one  Li 1  ,  two  I'  and  1.55  waters  in  a  cubic  simulation  eel! 
with  edge  length  1.862  nm  and  periodically  replicated  in  the  (x,y) 
plane.  The  metal  surface  at  z  —  0.932  nm  has  qM  =  +  |c| 
(anodic).  Figure  1  shows  the  density  profiles  averaged  over  the 
xy  plane  as  a  function  of  z  for  all  components  of  the  system,  flic 
iodide  distribution  is  sharply  peaked  near  0.7  nm  just  inside  the 
repulsive  portion  (dash  line)  of  the  9-3  potential.  The  iodide  dis¬ 
tribution  is  compact  and  the  Li 4  distribution  is  diffuse.  At  nega¬ 
tive  7.  the  Li  density  has  the  gradual  fall  off  expected  for  an 
electrostatically  bound  species.  The  main  component  of  the  Li 
density  ends  at  0.05  nm  which  we  interpret  as  the  position  of  the 
OUP,  being  approximately  two  water  diameters  from  the  repulsive 
wall  at  7  -  0.682  nm.  In  the  region  between  0.05  and  0.3  nm  the 
cation  has  finite  probability  of  penetrating  the  inner  layer.  This 
process  is  not  part  of  usual  Stern-Gouy-Chapman  theory,  'flic 
structure  near  the  metal  surface  in  the  densities  for  water,  11,  and 
PC  (point  charge  of  the  st2  model)  arises  from  an  oriented  layer 
of  water  with  a  PC  pointing  at  the  metal.  This  an  important  result 
of  our  model  with  implications  for  modulation  spectroscopy. 

Application  to  Surface  Spectroscopy 

The  idea  of  altering  the  position  of  the  outer  plane  to  change  the 
field  across  the  inner  layer  has  been  proposed  many  times  ’>  2. 
This  effect  can  be  modeled  with  the  primativc  system  used  hcic. 
Based  on  the  total  charge  distribution  given  in  Figure  I  we  see 
that  the  iodide  and  oriented  surface  waters  determine  the  charge 


density  near  the  metal.  Note  the  positive  charge  peak  at  0.65  nm. 
It  and  weaker  positive  peaks  near  0.4  and  0.2  nm  make  it  <IilTlctilt 
for  positive  ions  to  approach  the  interface  even  though  the  latter 
carries  net  negative  charge  because  of  two  adsorbed  anions.  Two 
opposing  effects  operate  as  cation  size  is  increased,  and  their  bal¬ 
ance  determines  the  field  cfTcct.  First  the  OIIP  will  move  to  more 
negative  7.  due  to  larger  radius  and  electrostatic  repulsion  from  the 
positive  charge  peaks.  Second  at  larger  radii  the  cation  hydration 
shell  is  softer  making  it  easier  for  the  ion  to  penetrate  the  inner 
layer.  In  Figure  I  the  tail  between  ().()  and  0.3  nm  indicates  how 
difficult  this  is  for  Li  ion  with  its  strongly  bound  solvation  shell. 
More  calculations  exploring  this  effect  are  in  progress  and  will  be 
presented  at  the  meeting. 

Conclusion  and  Acknowledgement 

Molecular  dynamics  calculations  show  that  a  simple  model  based 
on  sound  chemical  ideas  reproduces  phenomena  familiar  from  ex¬ 
periments  on  the  electrochemical  interface.  Important  new  in¬ 
sights  arise  concerning  the  time  averaged  electric  Helds  across  the 
inner  layer  which  in  turn  deepen  our  understanding  of 
spectroscopic  probes. 
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Figure  1.  Density  profiles  for  two  I\  one  Li*  and  155  st2  water 
molecules  next  to  immersed  electrode.  Anodicaily  charged 
electrode  qM  =  +  |e|.  Iodide  distribution  is  compact  and  Li+  dif¬ 
fuse.  Metal  on  rhs,  dielectric  constraining  boundary  on  Ihs.  Dis¬ 
tance  across  the  cell  1.892  nm.  Repulsive  portion  of  wall,  dashed 
line,  at  |z|  =  0.682  nm.  Temperature  290K. 


